Theory of measurements of electrodynamic properties in anisotropic superconductors 
in tilted magnetic fields. Part I: flux flow and Campbell regimes. 
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The vortex dynamics of uniaxial anisotropic superconductors in magnetic fields applied with arbi- 
trary orientation is theoretically studied. Focus is on the model for electrical transport experiments 
in the linear regime. Relevant vortex parameters, like the viscous drag, the vortex mobility and pin- 
ning constant (with point pins), together with the flux flow and Campbell resistivities, are derived 
in tensor form, in the very different free flux flow and pinned Campbell regimes. The applicability 
to the various tensor quantities of the well-known scaling laws for the angular dependence on the 
field orientation is commented. Moreover, it is shown that the experiments do not generally yield 
the intrinsic values of the anisotropic viscosity and pinning constant. Explicit expressions relating 
measured and intrinsic quantities are given. 
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I. INTRODUCTION 



Many superconductors of wide interest and recent discovery, such as iron-based superconductors 1 , MgB 2 2 and 
cuprate superconductors, 3 have in common an intrinsic material anisotropy, essentially uniaxial, arising from their 
crystal structure. The material anisotropy has a profound impact, among the others, on the vortex dynamics and 
on the related pinning phenomena. Such properties have been much studied due to their importance both for 
unraveling the fundamental physics of the underlying superconductor and in view of technological applications. As an 
example, recently a great deal of effort is devoted to the artificial tailoring of pinning on YBa2Cu 3 07-,5 based coated 
conductors, 4 ' 5 mainly through the introduction of extended defects whose effect is to further enhance the anisotropic 
behaviour. 

I focus on the model for electrical transport measurements in the linear regime in the mixed state, since this class 
of measurements is largely used in the study of vortex dynamics. 

The interplay between the material anisotropy and the preferential direction introduced by the magnetic field B 
determines a non-straightforward relationship between the applied current density J and the corresponding electric 
field E. Indeed, by applying with a generic orientations J and B, the vortices move under the effect of the Lorentz 
force J x £?$o and induce (Faraday's law) an electric field E which is in general not parallel to J even in isotropic 
superconductors. In anisotropic superconductors, additionally, in general vortices do not move parallel to the Lorentz 
force, further reducing the possibility to obtain an electric field E parallel to the applied J. 

As a consequence, the measured quantities, such as flux flow and Campbell resistivities (and their vortex coun- 
terparts, the vortex viscosity and the pinning constant), depend on the angles between B, J and the anisotropy 
axes. Hence they require a tensor representation and become not straightforwardly related to the material intrinsic 
properties. 

Previous works addressed some aspects of the problem, such as: anisotropic flux flow in the pin-free d.c. regime; 6-11 
pinning in non linear regimes in tilted fields, studied in the perspective of magnetization measurements; 12 ' 13 two- 
dimensional anisotropic pinning with isotropic viscous drag and fixed magnetic field orientation; 14 coupling between 
anisotropic two-fluid currents and vortex motion, the latter described within an isotropic framework. 15 ' 16 

In this work I propose a generalized treatment, centered on the force equation for the vortex motion, referring to 
uniaxial anisotropic superconductors in a magnetic field applied with generic orientation. Both the material anisotropy 
and pinning, the latter limited to weak random point pins only, are considered and studied in the very different free 
flux flow and pinned Campbell regimes, for arbitrary angles between B, J and the anisotropy axis. 

The objective is to provide all the quantities of interest as tensors, clearly distinguishing between intrinsic and 
measurable properties. Moreover, in several examples they will be cast into expressions directly exploitable in the 
experiments and applied to analyze experimental data. 

This work is organized as follows: in Section II, the electrodynamics model for the flux flow regime is recalled 
and recast in a suitable form; in Section III the vortex viscosity, vortex mobility and flux flow resistivity tensors are 
computed starting from the vortex force equation; in Section V the treatment is extended to the a.c. Campbell regime, 
yielding the various pinning tensors. Sections IV and VI are devoted to experimental aspects, providing examples for 
the measurement of the mixed state anisotropic resistivity and of the pinning constant, respectively. 



II. VORTEX MOTION ELECTRODYNAMICS MODEL 

In their works, 6,7 Hao, Hu and Ting (HHT) addressed the problem concerning the d.c. flux flow resistivity in 
anisotropic superconductors in the mixed state. The model holds in the linear regime in a homogeneous superconductor 
with an uniform magnetic field applied along a generic direction, in the London limit. The model neglects pinning 
and assumes vortices to be straight and rigid flux lines moving in a uniform current field density. This basic model 
does not take into account more complex phenomena such as helical instabilities of the vortex lines, 17 which can 
occur when there are current components parallel to the applied field, or flux-line cutting effects, 18 or the breaking 
of vortex lines into pancakes, occurring in the extreme anisotropy, layered superconductors. 17 Accordingly, I do not 
consider any electric field component || B. 12 An intrinsic flux flow conductivity tensor 19 a^g is introduced, 6 ' 7 relating 

the electric field E induced by vortices moving with velocity v, which can be expressed through the Faraday's law 
as: 20 

E = Bxv (1) 
to the current Jt which determines the flux flow dissipation, so that, within the linear response theory: 

J T = d { $E (2) 
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The inverse of crj^ yields the intrinsic flux flow resistivity tensor p K g : 



W. 



(3) 



It is important to note that and p^g are intrinsic quantities which in general arc not directly measured. Indeed, in 

typical experimental setups the imposed current J needs not to coincide with Jt of Eq. (3), so that the experimentally 
measured resistivity tensor pg, defined as: 



E = p ff J 



(4) 



will be in general different from p-g . The above equation states that in general, because of the anisotropy, E J. 

An important outcome of this paper will be to obtain explicit relations between the measured {pg) and the intrinsic 
{p^ff) nux fl° w resistivity tensors. 

The fact J ^ Jt can be understood considering that the electric field E induced by the vortex motion is by definition 
perpendicular to B 7 i.e. E ■ B = 0, and therefore: 6 ' 7 



J T ■ B = 



(5) 



The latter expression imposes that Jt cannot be freely oriented with respect to B, so that Jt is in general distinct 
from J, which instead can be externally applied with arbitrary orientation. One can write: 6,7 



J = J T + Js 



(6) 



where Js is a supercurrent density, parallel to B (i.e. J$ = JsB, being B the B unit vector), giving rise to no 
dissipation related to vortex motion (flux line cutting 18 is neglected). 

The explicit expression for the a^} has been previously derived. 8-11 I follow Ref. 8 which, working within the 
Time-Dependent Ginzburg-Landau theory (B < B C 2), accounts for both ohmic losses and order parameter relaxation. 
I first specify the frame of reference. The crystallographic axes are taken as the coordinate axes of a cartesian frame 
of reference (Fig. 1), so that x = a, y = b and z = c, being the latter the axis of the uniaxial anisotropy. In this frame 
of reference the phenomenological mass tensor, 21 ' 22 which can be used to describe in the London limit the material 
anisotropy, is diagonal: 






' 1 







= m ab 
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• 




Vo 








mM 



(7) 



having defined the in-plane mass m a b = m, the out-of-plane mass m c and the anisotropy factor 7 2 = m c /m a b. 

The mass tensor M contains all the information concerning the material anisotropy, since the only source of 
anisotropy that will be considered is the effective mass of the charge carriers (this implies, for example, that the 
possible anisotropy of the scattering time of the normal carriers is neglected). By neglecting the Hall contribution 




FIG. 1. Principal frame of reference. The magnetic induction field B is also depicted, applied along a generic direction at the 
polar 9 and azimuthal cj> angles. 



and assuming the same anisotropy axes for the normal resistivity and mass tensors p n — p n ,uM (i.e. the mass tensor 
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is the only source of anisotropy also for the normal state), the tensors p~g and (jg 1 are also diagonal in the stated 
frame of reference. Therefore one has: 

/pguGM) o o \ 

pf{B,6)=\ pf tll (B,0) (8) 
V P { £ 33 (B,9)J 

Neglecting the weak field dependence of p n and writing down explicitly all the field dependencies, one obtains: 6,8 

P%{B,6)/p n ^=F{B/B c2 {6)) (9) 

This equation implies two relevant results. First, it is consistent with the important scaling law 22-25 according to 
which the field dependence of the anisotropic physical quantities in the mixed state depends on the B/B c2 {9) ratio 
only. Second, it shows that all the three tensor elements of p^ share the same field dependence through a common 
function T(BjB c2 ). 

Therefore Eq. (8) can be rewritten as: 

p ( ff\B,6) = PnA1 T(B/B c2 (6))M = 

= pf^{B/B c2 {6))M (10) 

The above Equation highlights the important property that the field and angular dependence of the whole intrinsic flux 
flow resistivity tensor can be represented by a single scalar function, namely the element Pff > 11 (B/B c2 (9)). Similarly, 
the conductivity tensor is: 

3$(B,e)=a£ 11 (B/B c2 (e))M- 1 (11) 

In the following, for ease of notation the angular and field magnitude dependence will be explicitly written only in 
the Equations reporting the main results. 



III. FLUX FLOW REGIME 



In this section the well-known vortex force Equation, 26 ' 2 involving the balance of forces acting on an individual 
vortex, is studied in the regime of pure flux flow. This regime is realized by an ideal pin-free superconductor, 
indifferently in d.c. and a.c. regimes, whereas in real superconductors it can be obtained with d.c. currents sufficiently 
large to overcome the pinning forces, or with a.c. currents at high enough frequencies. This treatment will provide a 
generalized formulation which allows (i) the determination of the vortex viscosity tensor fj (Sec. Ill A), which gives 
many physical insights, and the related vortex mobility t2 v tensor (Sec. IIIB); (ii) the subsequent computation of the 
measurable flux flow resistivity tensor pg (Sec. IIIC); (hi) the basis for the extension to a.c. regimes with pinning, 
such as the Campbell regime treated in Sec. V. 

In pure, stationary, flux flow regime, with no pinning and hence no creep contributions, the viscous drag force (per 
unit length) fjv acting on the vortices moving with velocity v exactly balances the Lorentz force Fl = &oJ x B so 
that the force Equation, written in vector form, is: 

t]v = ^qJxB (12) 

In the above Equation, it is apparent that Fl and hence the viscous drag tensor fjv belong to the plane perpendicular 
to the magnetic field. 



A. Vortex viscosity tensor 

To obtain fj using Eq. (12), I consider an arbitrary field B = B(sm8 cos <fi, sin6> sin0, cos#) (see Fig. 1) and a 
generic current density J. Vortices move under the Lorentz force <&oJ x B. Using Eq. (6), J x B = Jt x B since 
Js || B. Then, Jt can be expressed in terms of v through Eqs. (1) and (2), so that, using Eq. (11): 

f?=*o(ffj ) (Bxil)) xB = -^oBB x B i g ) B x v = 

= -$ Q BB x a i £ 11 M- 1 B x v = ^ Q Ba^] u Mv (13) 
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valid W. 

In the above, the tensor B y comes from a general property of the cross product, 28 according to which, given the 
generic vector a, one has B x a — B x a with B y defined from B as described in Appendix A 1. Moreover, for the sake 
of compactness, the tensor M has been introduced: 



cos cos V sin ( 



M(0, 



ByM~ By = 



'cos 2 9 + 7 2 sin 2 <f> sin 2 9 -7 

2 cos sin </> sin 2 9 cos 2 + 7~ 2 cos 2 0sin 2 6* — sin cos 6* sin 9 



"7" 

— cos </) cos 6* sin 



• cos sin <p sm W 

-2 „„ c 2 ^„:„2 , 



(14) 



— sin 6 cos 9 sin 



sm 



By equating the first and last member of Eq. (13) Vw, one obtains: 



fj(B,9,0')=<f Q BB x B { ff ) B x = 

= <S> Ba { j ; ) 11 (B/B c2 (9))M(9,ct>) (15) 

It can be noted that the viscosity tensor fj does not obey to the angular scaling law and that it depends also on 0, 
even if the superconductor is uniaxially anisotropic, because of the Faraday-Lorentz contribution included in M(9, (j)). 
Analogously to and an "intrinsic" viscosity tensor (diagonal like a^) can be introduced as follows: 

r} W (B,9) = <S>oBa$(B/B c2 (6)) = 

= V ^(B/B c2 (9))M- 1 (16) 

which obeys the angular scaling law and whose elements are r$ — ^oBa^ u . Equation (15) can be thus rewritten 
as: 

fj(B, 9, <j>) = -S x (9, <p)f}^ (B/B c2 {9))B x (9, 0) 

= rj^(B/B c2 (9))M(9^) (17) 

Before computing pjj, some relevant properties of the viscosity tensor will be highlighted. 



1. Eigenvectors and rank of the viscosity tensor 



The eigenvectors of the viscosity tensor fj are interesting since, by definition, they are specific vortex velocities v e 
for which the viscous drag force is parallel to v e , f\v e = i] e v e , where the scalar r\ e is the eigenvalue of v e . Hence, for 
these special directions the viscosity tensor behaves as a scalar, similarly to the behaviour in isotropic materials. 

The eigenvectors and eigenvalues can be computed through standard tensor algebra, 29 yielding the results reported 
in Table I. 

In the Table the function e 2 {9) has been introduced: 

e 2 {9) = cos 2 (9 + 7~ 2 sin 2 6> (18) 

which corresponds to the well-known angular-dependent anisotropy parameter 25 ' 30 e 2 (9), which defines, among the 
others, the angle dependence of the critical field B c2 {9) = B c2 (0)/e(9). It can be seen that the trivial v e s || B has 



Eigenvector v e ,i 




[J,Ve] 

for 4> — 


7T 

2 


l:(sin0, — cos</>, 0) = (B x z)/\B x z\ 


(«) 2 

mi* 


[B x x, 


X] 


2:(cos 9 cos 4>, cos 6 sin <f>, —sin8) = Bxv s i 


ml 


[x, B x 


X] 


3:(sin#cos0, sin # sin 0, cos 6<) = B 





N.A. 





TABLE I. First two columns: eigenvectors and corresponding eigenvalues of the viscosity tensor within uniaxial anisotropy; 
third column: eigenvectors computed for (j> = tt/2 and corresponding current direction forming a right-handed orthogonal basis 
with B and v e (see text). 
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zero eigenvalue (as expected from expression (14)) and that the other two are both _L B. Hence the three (unit) 
eigenvectors of Table I constitute an orthonormal right-handed basis "d" , i.e. they are the principal axes for the 
viscosity tensor. The viscosity tensor fj is then by definition diagonal in the frame of reference with basis "d" . 

The transformation matrix for this new frame of reference is D = \v e i v e2 t^], 29 i.e. its columns are given by 
the coordinates (in the standard frame of Fig. 1) of the three unit eigenvectors. In the following, the quantities 
represented in this frame of reference will be denoted with the superscript "(d)": for example, given a generic vector 
a and matrix A, a* d ) = D^a and A^ = D~ 1 AD. 29 One obtains: 

(e\6) (A 
fj( d \B,6)=r ] { ?}{B/B c2 {6))\ 10 = 

\ 0/ 

= f 1 ^(B/B c2 {6))M {d \e) (19) 

It can be seen that fj^' is a symmetric non diagonal matrix, with zero elements in the third row and third column. 
The zeroes in the third row ensures that the viscous drag force component || B is zero, as anticipated in the comment 
of Eq. (12). The zeroes in the third column makes the viscous drag force independent from a possible component 
of v parallel to B. Thus, the force Eq. (12) with the results (15) and (19) imply that v ± B, consistently with the 
framework of the present model in which a motion of rigid vortices parallel to their axes is meaningless. 

Moreover, it is evident that the viscosity tensor, defined within the three-dimensional space, is not invertible, since 
its rank=2. Hence the viscosity tensor is a bijection not in the whole three-dimensional space but only within the 
plane perpendicular to B. This implies that, for a given Lorentz force, the corresponding vortex velocity is univocally 
determined, property which is useful in the computation of the vortex mobility tensor. 

Going back to the eigenvectors, since v e \ and v e2 are _L B, the force Equation (12) implies that for both of them it 
is possible to find a proper current direction J which gives a right-handed orthogonal basis [J, B,v e ], which evidently 
ensures that E || J. 

Taking B in the y-z plane (</> = tt/2) for the ease of notation and without loss of generality, these current directions 
are reported in the third column of Table I and represented in Fig. 2, where the left and right panels corresponds 
to the second line and first line in the Table, respectively. In the Figure, yd — (B x x) \\ J has been introduced 
for later reference. It is worth stressing again that in general E J, which makes more complex the experimental 



z z 




FIG. 2. Specific field-current configurations (see text). 

determination of the material resistivity. In this view, the above two special current-field configurations are especially 
useful. 



2. Vortex viscosity along the principal axes 

I now consider specific current-field configurations where both vectors are directed along the principal axes of a 
uniaxial superconductor. I compare the viscosity obtained from the tensor analysis to the results 11 computed on the 
bases of the Bardeen-Stephen 31 model. From Table I and Eq. (Ill A) , it is easy to check that for these configurations 
one has the usual right-handed orthogonal basis [J, B,v] and that f]v — r\ e v. As it can be seen in Table II, the results 

obtained through the two above mentioned procedures are coherent. In the Table the relation ^^(f) = (0)7, 
easily obtainable within the Bardeen-Stephen model, is used. 
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Set 
[J,B,v] 


angles 


n from Ref. 11 


r? e from Table I 
and Eq. (Ill A) 




x,z,y 
x,y,z 
z,y,x 
z,x,y 


0,.. 

7T 7T 

2 ' 2 

7T 7T 

2 ' 2 

- 

2 ' u 


( c ) 

(6) (a) ^ (c) 

(f>) (c) -1 
Va ~% 7 
(a) (c) _l 

V b ~V b 7 


^(0) 

^ ) (f)=^i ) (0)7 
4 ) (f)=^l ) (0)7- 1 

^ ) (f)=^i ) (o) 7 - 1 


diag(r,«(0) 1?? «(0),0) 
diag(^(f), 0,^(1)) 
diag^gd), 0,^(1)) 
diag(0, ^(f), ^(f)) 



TABLE II. First two columns: geometric configuration; third column: n from Ref. 11, following the notation used there (i.e. 
subscript and superscript denote velocity and field directions, respectively); fourth column: the scalar viscosity n e obtained 
from Table I and Eq. (Ill A); fifth column: full fj tensor. 



One additional comment can be done by considering the full fj tensor explicitly computed for the field orientations 
along the principal axes and reported in the fifth column of Table II. In this case, the viscosity tensor is diagonal: 
this means that, by taking B and J along two principal axes, the vortex motion (both in terms of v and Fl) occurs 
along the third principal axis. 

Moreover, when B is parallel to a principal axis, it can be easily verified that the viscosity tensor fj is diagonal if 

and only if a^g is itself diagonal, meaning that the same principal axes of are inherited by fj. This property will 
be exploited when treating the pinning constant tensor in Sec. V. 



B. Vortex mobility tensor 

The force Equation (12) can be written in terms of the vortex mobility tensor p v as: 

v = Jl v <5>oJxB (20) 

Since the vortex viscosity tensor has rank=2, i.e. it is not invertible, p v cannot straightforwardly be derived as the 
inverse of fj, as the relation between the corresponding scalar quantities would suggest. Thus, I exploit the results of 
Section III A: within the plane _L B, ji v is the inverse bijection of fj whereas the vortex mobility for the (physically 
not existing) Lorentz force component || B can be safely set = 0. Hence, in the frame of reference "(d)" the vortex 
mobility is the following diagonal matrix: 

(^'diag^ 2 ,!^) (21) 

In the principal frame of reference, one obtains: 

fiv = Dfl^D- 1 (22) 



C. Measured flux flow resistivity tensor 

The flux flow resistivity tensor pg can be now computed in terms of the vortex mobility Ji v . By left cross- multiplying 
both members of Eq. (20) by B one obtains, after a little algebra, E = —B x fl v $ BB x J, valid V J, which implies: 

p ff = -B x fi v <$> BB x (23) 

It is of paramount importance that Eq. (23) represents the measured, apparent, flux flow resistivity tensor, different 
from the intrinsic flux flow resistivity tensor. The relation between pg (measured) and p^g (intrinsic) can be made 
explicit as follows. In the frame of reference "d", considering that D~ 1 B X D = z x and using Eq. (21), one obtains: 

p ff = ^ Bb(-i x ji^i x )b- 1 = 

=<M?i> diag( v: 2 '%- i = 

, Dfj^D^ 1 <5> Bfj 
= ®oB 



»7ii e 



2 



2 
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Using Eqs. (16) and (III A) the role of a^j) emerges: 



so that, in terms of p^ — ■ 



MB,0,j>) =p$ A1 (B/B c2 (6))^^ (25) 

It can be noted that pg, like the viscosity tensor fj of Eq. (15), does not obey to the angular scaling law. This 
expression for the measurable flux flow resistivity tensor is particularly useful and as such will be largely used and 
extended in the rest of this work. It can be checked that, as shown in Appendix B, the same result can be obtained 
following the whole HHT model but, as already mentioned, the present derivation from the force Equation makes 
possible extensions to a.c. regimes as it will be done in Sec. V. 

It is worthwhile to stress that Eq. (25) implies that the measured flux flow resistivity tensor pg is symmetric 
(respecting the Onsager principle 32 ) but non-diagonal even if one has chosen to neglect the Hall terms in the intrinsic 
(diagonal) tensor p~g . This important property describes the physical phenomenon that, by applying with generic 
orientations J and B to a superconductor (either isotropic or anisotropic) , the vortex motion induces an electric field 
E which is in general not parallel to J. Moreover, in anisotropic superconductors vortices do not move in general 
parallel to the Lorentz force direction J x B, further reducing the possibility to obtain E parallel to J. 



IV. DETERMINATION OF THE INTRINSIC RESISTIVITY TENSOR ELEMENTS FROM 

EXPERIMENTAL MEASURED QUANTITIES 

In this section I relate the experimentally obtainable quantities to the intrinsic resistivity tensor, in various typically 
accessible experimental configurations. The results, Eqs. (27)-(32), give the tools to extract the tensor elements Pff,u 
from the experiments. The issue of the true current path is of paramount importance, and it has to be solved 
separately. 33 Here, I focus on the scalar resistivity p^ measured along the direction of the applied current J, which 
is computed in general as: 



J-E/J 2 = J-E/J = (^pfj-J (26) 



where in the flux flow regime p = pg (measurable). 

Since E _L B, then J ■ E — Jt ■ E 1 coherently with the fact that Jt is involved in the vortex motion dissipation 
(Eq. (2)). In the above equation, J • E is the electric field component, parallel to J, which has to be experimentally 
measured in order to determine p^ . 

First I consider a static field B = B (sin 9 cos (j), sin 9 sin (f>, cos 9) with generic orientation (depicted in Fig. 1) and 
an external current J = Jx applied along the x-axis. This configuration is commonly employed to obtain "a-6 plane" 
properties. 

The resistivity p^ x \ where the superscript "(a;)" refers to the orientation of the applied current, can be computed 
through Eq. (26) with p — pg given by Eq. (25), yielding: 

m Pffii (°) sin 2 9sin 2 (f) + p^\J9)cos 2 9 

p^(9,ep)=p { j F \j9 y ff ' 1 j; ' r ; ff ' 331 ' (27) 

Pff,ll^> Sm U + Pff,33VV) COS V 

This result, here obtained through tensor algebra, coincides with Ref. 6 and 7. Equation (27) shows the very 
important result: the measured resistivity is in general a mixture of the intrinsic flux flow resistivity tensor elements 

(in-plane or a-b plane resistivity) and Pg\ 3 (out-of-planc or c-axis resistivity), even if the external current flows 
parallel to the a-b planes. 

The intrinsic resistivity p'gn (9) can be directly accessed by making measurements within the well-known "maximum 
Lorentz force", where J = Jx as before and in addition B(6,(j) = ir/2) = _B(0,sin#,cos#) g y-z plane, so that the 
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magnetic field forms, as demonstrated in Sec. IIIA1, a right-handed orthogonal basis [J = x, B,v = Ft] with the 
current and the vortex velocity (the latter being an eigenvector of fj). 

This configuration is represented in the left panel of Fig. 2. One obtains: 

pU(e,w/2) = p%] 11 (0) (28) 

This is the only easily accessible configuration where one has direct access to the a-b plane intrinsic flow resistivity. 
Another typical configuration is J || z, 34 often used to study Pff] 33 - The result is, independently from the angle 4>: 

p%\ A6)p ( l\j 6) sm 2 6 

p^(6) = ... ff * ,IJ ff^) (29) 

^(6)^6 + ^(6)^6 

which yields Pff\ 3 only for p( z \-n/2) = p 33 (tt/2), i.e. for B G x-y plane. Observing this Equation, it can be seen that 
it is not possible to directly measure Pff\ 3 (6) for each 6, contrary to what happens for p^\ 1 (6). 

Experimental results analogous to those represented by Eq. (29) can be obtained by taking J = Jx and B(8, (f> = 
0) = B(smd,0,cos8) G x-z plane, so that the magnetic field forms the (varying) angle 6 with the current J: 

P^\A6)p ( Jf\,(6) cos 2 8 

p& (6, 0) = 0%33V_j (3Q) 

^(6)^8 + ^(8)^8 

This is another common experimental configuration, 35 which again yields an experimental resistivity arising from the 

admixture of p^gn and p^g 33 - 

Indeed, with a little algebra (see Appendix C) it can be demonstrated that the direct measurement of p^g 33 (6) 

is impossible in general: whichever combination of directions of B and J is chosen, pff\ 3 (0) cannot be directly 
accessed: the measured angular dependence of the vortex-state, c-axis resistivity in tilted fields never corresponds to 
the material-dependent quantity. 

Hence, p^ 33 (6) must be indirectly determined through coupled, complementary measurements. An obvious choice 
for one of the two needed configurations is the maximum Lorentz force configuration (Eq. (28)), while the other could 
be equivalently chosen between the two last commented configurations (Eqs. (29) and (30)). Another possible choice 
is the configuration discussed in the study of the viscosity tensor eigenvectors (see Sec. Ill A 1): B(6, ir/2) and J G y-z 
plane along the yd = B x x direction, as depicted in the right panel of Fig. 2, yielding the right-handed orthogonal 
basis [J = yd,B,v = Fl =x]. Although it is seldom used in the experiments (requiring samples grown with the c-axis 
parallel to the sample surface 36 ), it will be useful in the subsequent discussion about the Campbell resistivity tensor 
(see Section V) since it ensures E || J contrary to the two configurations of Eqs. (29) and (30). 

In this case one obtains: 

p(«-)(g) = „ P^ e )Pffff) (31) 

P { ff] n (0) sin 2 6 + p^ 33 (6) cos 2 6 

By combining Eqs. (28) and (31), the intrinsic p^ 33 (6) can be written down in terms of the measured resistivities: 

(0 ta ^_ P {yd] '(g) P {x) (A) sin 2 6 
PffM") - p {v d ) ( ) COS 2 W 

An analogous result can be obtained by using Eq. (29) or (30) instead of Eq. (31). It is also worth noting that 
the above Eq. (32), because of the subtraction in the denominator, 37 yield an indeterminate form 0/0 for 6 = 0, 
implying that p^ 33 (8 = 0) cannot be measured. This is true in general: it can be easily seen by writing down the 

measurable resistivity tensor (Eq. (25)) for 6 = 0, obtaining pg(8 = 0) = diag(p^' ) 11 (0), p~g xl (0), 0). 38 Since p^ 33 (0) 
does not appear in it, no current direction can yield the c-axis resistivity for 8 = 0. Hence the material resistivity 
along the anisotropy axis, i.e. the intrinsic c-axis resistivity, cannot be measured if the magnetic field is parallel to 
the same axis, not even by resorting to the measurement of the transverse resistivity components (i.e. through the 
determination of the transverse components of the electric field). 

On the other hand, for a generic ^ it has been shown that the correct process to determine the intrinsic c-axis 
resistivity p^ 33 (9) necessarily requires the combined measurements of the resistivities with two distinct current-field 
setups. 
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V. CAMPBELL REGIME 



A totally different limit with respect to flux-flow is given by the pinning regime. There, the d.c. response is zero. 
However, the a.c. resistivity is well measurable (Ref. 27 and references therein). To my knowledge, no complete 
treatise of the anisotropic resistivity in the pinning regime exists. In this Section I consider the vortex motion 
under a.c. currents flowing in anisotropic superconductors with random point pinning centers, described in the weak 
collective pinning regime. 25 In particular, I focus on the well-known Campbell regime, in which the pinning action on 
vortices dominates over the dissipative viscous drag and the pinning force is proportional to the small displacement 
u of the vortices from the pinning centers. The above two conditions are achieved respectively for frequencies smaller 
than the so-called pinning frequency (of the order of a few MHz in conventional superconductors 26 and of several GHz 
in high-r c superconductors, 27,39-41 low-T c thin films 42 and more complex superconducting heterostructures 43 ) and for 
currents small enough to ensure the validity of the linear regime. Considering only point pins (i.e. zero-dimensional 
pinning centers) ensures that no further preferential directions are introduced in the superconducting system, in 
contrast with extended (e.g. linear or planar) defects. 

The starting point is the force equation written in the sinusoidal regime e 1 "*, including the pinning force and 
neglecting the viscous drag: 

$ J x B (33) 

MJJ 

where k p is the pinning constant (also called Labusch parameter) tensor and u = v/(ico) is the displacement of the 
fluxon from the pinning center. No creep phenomena are considered. 

In the Campbell regime, losses due to finite real conductivity are neglected, so that the resistivity is purely imaginary 
and defined, in isotropic superconductors, as: 

Pc = —r—u = w MoA c (34) 

where Ac is the Campbell penetration depth. 44 

Going back to anisotropic superconductors, it is evident that the force Equation (33) in the Campbell regime is 
formally equivalent to the force Equation (12) written for the pure flux flow regime. Hence, by extending to the a.c. 
regime the electrodynamics model for the electrical transport in the mixed state extensively discussed and exploited 
in Sec. Ill, it is straightforward to introduce the series of tensors [k p /(iu), k p ^/{\uj), — itr^, ipg* , ipc] which are dual, 
both in terms of roles of the tensors and of their relationships, to the already studied series [fj, fj^\ p-g , pg ]. 
Therefore the following expressions, analogous to Eqs. (Ill A), (16), (3) and (B4), hold: 



kp 


= -bj^b x 




= uj$ Ba 


(i) 
C 






1 


PC 


= bA- 


=(») 

Pc 





{pfB) 



B, 



■ B 



(35a) 
(35b) 

(35c) 
(35d) 



where the latter expression holds for a diagonal p^\ as will be taken in the following. It is worth noting that the above 
Equations do not complete the model in the Campbell regime, since an explicit expression for the intrinsic tensor 
Pq 1 is still needed. Whereas in the flux flow regime the problem about the corresponding intrinsic tensor p^ was 

addressed by the TDGL treatment, 8 the determination of the intrinsic tensor p^) will be addressed in the following 
Section. 



A. The Campbell resistivity tensor elements 

I assume that p^) is diagonal like p^}: this ensures that, according to the final comment of Sec. Ill A 2, choosing B 
and J along two principal axes, the corresponding pinning force (and vortex velocity) will be parallel to the third axis, 
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as it is reasonable to expect. This assumption does not define the whole tensor, since it does not give any information 
about the two (within the uniaxial anisotropy) non-zero diagonal elements, Pq\i{@) and Pc\ 3 (9). 

As already seen for its flux flow counterpart, Pc\ 3 (9) is not directly accessible through real current- field configu- 
rations. Therefore, it has to be derived indirectly starting from two experimentally measurable resistivities. I choose 
P < q\6) and p i Q d \9), which are defined within the two specific field-current configurations (<fi = n/2) commented in 
Sec. IV and depicted in Fig. 2. It is worth recalling that, in the present Campbell regime, the applied J is a 
(low- frequency and small intensity) a.c. current. Rewriting Eqs. (28) and (31) for pc, one has: 

p£\ 1 (B,9)=pP(B,8) (36a) 
Tg-v ' - cos 9 

According to Sec. Ill Al, both the field-current configurations "(x)" and "(yd)", related to p^ and p^ d \ correspond 
to eigenvectors of the vortex velocity v, and therefore of the vortex displacement u \\ v. Hence the pinning force 
<e u is || u and is completely described by the scalar eigenvalue pinning constant k p . e , which within the 

(x) (v<i) 

two configurations will be denoted as k p or k p . This property explains the present choice of the field-current 
configurations, since other geometries would have given different orientations for u yielding k p u u and thus preventing 
a description of the pinning force through one scalar k Pye only. 
The definition (34) allows then to write down the following: 

#™-iF%ir <37a) 

The actual values of the pinning constants are determined as follows. 



B. The pinning constant tensor elements 

By simple physical arguments, 45 for an isotropic superconductor one can evaluate the pinning constant k p by 
equating the maximum pinning force k p r p acting on a vortex, where r pin denotes the action range of the pinning 
centers, with the maximum Lorentz force J c $o exerted when the current equals the critical current density. Since 
for core-pinning r p i„ ~ £, i.e. the coherence length which define the length scale for the spatial variations of order 
parameter and therefore the radius of the vortex core, one can write 

k P = (38) 

where c ~ 1. 

For anisotropic superconductors, Eq. (38) must be specialized for the two current configurations "(x)" and "(yd)" 
used in Eqs. (37). The pinning constant k p x ^ [k p Vd ^j describes the pinning action on vortices moving along the ya 
[x] direction under the action of a current || x [ || yd}- hence the critical current density || x [Jc Vd ^ || yd], and 
r P in ~ ^ Vd ^ [ r P in ~ i-e. the coherence length along the direction y d [x] of the vortex movement. Therefore: 

-^^S^r (39a) 
^ HBJ)= *c4^«) (39b) 

One should note that fcp yd \ despite being referred to a current-field configuration of difficult realization in the exper- 
iments, is related to quantities {ji Vd ^ and t;^) which can be expressed within available theories. In this sense, Eq. 
(39b) gives a tool for subsequent elaborations. 
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Indeed, the above equations can be further developed using the Blattcr-Geshkenbein-Larkin (BGL) scaling law: 23-25 
in the London approximation, a thermodynamic or intrinsic transport property q of a uniaxially anisotropic supercon- 
ductor depends on the applied field B and angle 9 only through the ratio B/B C 2(9), so that q(B,9) — s q q lso (Be(9)), 
where q tso {B) is the field-dependent quantity in the equivalent (according to the scaling rules) isotropic superconduc- 
tor, and s q is a scaling factor typically equal to 1, 7 or e ±1 (9). 

The current densities and Jc Vd ^ are provided in Ref. 25 for point-pinning in the single vortex and small- 
bundle pinning regime (the scaling theory does not describe the large bundle pinning regime, in the highest field and 
temperature regions). It is found 25 that they scale as: 

4 X \B,6) = J l c so (Be(9)) (40a) 
J^\B,6)=e(6)Jr(Be(9)) (40b) 

where the full expression of Jl so (B) is reported in Ref. 25 for different pinning regimes. 

It is worth stressing that the critical current ,P c so {B) (and the corresponding kJf°(B) = c§oJ* so (B)/!; lso ) for the 
equivalent isotropic superconductor is not angle-dependent because point pins only are considered, whereas extended 
pinning centers would introduce preferential directions and therefore angle-dependent quantities even in the isotropic 
superconductor. 

Geometrically, the quantities £}- Vd \9) and £}- x \9) represent the maximum distances between the axis of the displaced 
vortex and the point pin. These distances are univocally defined since in the present model vortices are straight lines 
and the pinning centers are points with zero spatial extension. They scale as: 25 

f (w) (0)=e(0)f (41a) 
&\0)=i (41b) 

where £ = £ lso denotes the in-plane coherence length. 
By using Eqs. (40) and (41) in Eq. (39), one can write: 

k^(B,9) = s kix) k; s °(Be(9)), = l/e{6) (42a) 
k^(B,6) = s kiyd) k; s °(Be(9)), = e{9) (42b) 

It is worth stressing that, thanks to Eq. (36a), it is 

k$ 1 {B,0) = kW{B,0) (43) 

hence the in-plane pinning constant kp^ u (0) obeys the scaling law given by Eq. (42a). 



C. The Campbell regime tensors 

According to Eqs. (37a) and (42a), one obtains: 

p% u (B,6) = s PCtl JS°{Be{9)), 8poA1 = 1 (44) 
On the other hand, p^ 33 , using Eqs. (36), (37) and (42), becomes: 

(i) pgnCg^sin 2 ^ _ ^(£,0) sin 2 _ 

P CM"> > k^\ m _ cQs29 - 6 2 (0)-COS 2 



= p% n (B/B c2 (6)h 2 (45) 
Therefore, similarly to of Eq. (10), one can write: 

p<$(B,6)=p% n (B/B c2 (6))M (46) 

This is an important result of this paper: for point pinning, the anisotropy of the pinning tensors and the flux flow 
tensors is the same and is completely described by the mass anisotropy tensor. This result is very important also in the 
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study of regimes where both losses and pinning phenomena are equally relevant, placing themselves in a intermediate 
situation between purely dissipative flux flow and purely reactive Campbell regimes. These regimes, occurring at 
higher frequencies including the microwave range, will be studied in a future work. 46 

Lastly, due to Eq. (45), the intrinsic pinning constant tensor also satisfies the scaling law: 

k^(B,9) = k p %(B/B c2 {9))M^ (47) 
while the measured pinning constant tensor does not: 

l p (B,6) = -B X %W§ X = B c2 (6))M{8, 4>) (48) 

because of the additional contribution given by A4(0, <j>) (Eq. (14)). 



VI. APPLICATION TO EXPERIMENT: A.C. LINEAR SUSCEPTIBILITY MEASUREMENTS 

ANALYSIS 

In this Section I provide an example showing the additional information which can be gained using the results of 
the present work by performing an analysis of experimental data. Various techniques can be used to explore the 
Campbell regime, including inductance measurements, 47 a.c. linear susceptibility measurements 48 and the vibrating 
reed technique. 49 The pinning constant can be also determined through microwave measurements, 41 ' 50 but they 
typically go beyond the pure Campbell regime and, as such, they are beyond the scope of this work and postponed 
to a future study. 46 In the following, I focus on a.c. linear susceptibility measurements on YBa2Cu307_,5 samples in 
the mixed state, performed by varying both the direction and the intensity of the applied d.c. magnetic field. The 
source of the experimental data is Rcf. 51, where full details about the experiment can be found. The measured 
squared real part of the penetration depth A^j = X\ + X 2 C (Xl is the London penetration depth) is related to the 
squared Campbell penetration depth X 2 -, which, given the geometry of the experiment, is Xc.w oc k'^. The examined 
sample are two twinned YE^CusOy-^ single crystals, one irradiated at +30° with respect the c-axis in order to create 
columnar defects capable of reinforcing the pinning properties and one, virgin, used as reference. As a consequence, 
four sources of pinning are expected: point pins and three types of correlated pinning centers, namely the twins 
along the c-axis direction, the a-b planes and the columnar defects at +30° (in the irradiated sample only). It was 
concluded that pinning was stronger with the field aligned with columnar defects and, from a qualitative analysis, that 
an enhancement of pinning existed even far from the track direction. I show in the following that such findings can be 
put on solid quantitative grounds by exploiting the present model. In particular, I will use the model developed in the 
previous Section to remove the anisotropic response due to the material anisotropy and to point pins, which can give 
rise to a significant angle-dependent contribution. The latter can in principle mask the angle-dependent contributions 
arising from correlated defects. 

Here I consider the data taken from Fig. 5 of Ref. 51 for the irradiated sample only, where X 2 R (B) measured at 
various angles 9 at T=90.5 K is reported. In order to extract k Pt u(B,9), I extrapolate X R (B — > 0), take it as an 
evaluation of X 2 L and, by neglecting pair-breaking effects, estimate X 2 J (B) — X 2 R {B) — X 2 R {B — > 0). Then the pinning 
constant is obtained, using Eq. (34) with B = /iqH and normalizing with . The normalized pinning constant, 
kp,u/k™?i, is then free from geometrical factors. 

The result is reported in Fig. 3a. The pinning constant in the irradiated sample with B along the columnar 
defects (9 = +30°) is larger than along the other directions. Moreover, by examining the curves with 9 = —30° and 
9 = —70°, it was noted 51 that, despite being very similar at low fields, the two curves depart near half the matching 
field i?$ s=s 350 Oe. Thus, it was inferred 51 that the defects still determined an appreciable contribution to pinning 
even for the field tilted at 9 — —30°, 60° far from the tracks directions. Actually, this observation can be further 
substantiated by resorting to the result of Section V, namely Eqs. (42a) and (43): were only point pinning present, 
the in-plane pinning constant at different fields and angles should scale so that by plotting e(6)k Pi n vs Be{9) all 
the curves should overlap. Any possible residual difference between the scaled curves should then be ascribed to the 
effects of correlated defects only, since the scaling removes the background effect of point pinning centers. The result 
of such scaling is reported in Fig. 3b, where 7 = 7 has been used (taking 3 7 = 5-^8 does not change the result). 
It can be seen that, even by having removed the effect of point pins, along the tracks directions pinning remains the 
largest. On the other hand, the absolute values for 9 = —30° are now larger than those for 9 = —70°. Since at these 
angles the contribution of twins and a-b planes can be neglected (see also Fig. 2 in Ref. 51), one can infer that at 
9 = —30° the artificial defects still reinforce pinning, yielding a higher pinning constant than the one measured at 
9 = —70°. This result further substantiates, on quantitative grounds, the observation done on qualitative grounds, 
i.e. based on the field dependence, reported in the original work. 51 
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FIG. 3. Normalized pinning constant, as is (a) and scaled (b), as a function of d.c. applied field at selected angles, (a): bare 
quantities; (b): scaled quantities. Data digitized from Ref. 51 as reported in the text. 



VII. SUMMARY 



The electrodynamics model for transport measurements in the mixed state in uniaxial anisotropic superconductors 
has been discussed in the very different free flux flow and pinned (with point pins) Campbell regimes. Vortex 
parameters, like the viscous drag, the vortex mobility and pinning constant, have been derived in tensor form for 
arbitrary field orientations. It has been shown that the tensors describing point pinning share the same structure 
of the flux flow tensors, and that the measured quantities, differently from the corresponding intrinsic quantities, in 
general do not satisfy the angular scaling laws. It has also been shown that the experiments do not generally yield the 
intrinsic values of the anisotropic viscosity and pinning constant. Explicit expressions relating measured and intrinsic 
quantities have been given. An example of data analysis based on the results obtained in the Campbell regime has 
been provided. The results here obtained will prove necessary to investigate the angular dependence of the complex 
a.c. resistivity. 46 
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Appendix A: Tensor Algebra Notes 
1. Cross product 

The cross product a x b between a generic vector a and another vector b can be also computed as the application 
of a tensor a x to b: 

a x b — a x b 

where a* is defined starting from a : 28 





( o 


-a z 


a y 


a x := 


a z 





—a 




\-a y 


a x 






2. Factoring out tensors 

Here the following property is demonstrated: 



x (Ac) = A A- X (b x c) 



(Al) 



which holds V b, c and with A = diag(an, a 2 2, 033) an invertible diagonal tensor. 
Since it is easy to verify that 



(I?) x (AS) 



anx a 2 2U a 33 z 
a\\b x a 2 2b y a33b z 
a\\c x a 2 2C y a 33 c z 



= 011022033 



x y z 
b x by b z 

C X Cy C Z 



A 



(bxcj 



the thesis straightforwardly follows. 



Appendix B: Measurable flux flow resistivity tensor according to HHT model 

In Rcfs. 6 and 7, HHT provide an elegant and compact way to derive, from the intrinsic flux flow resistivity 
reviewed in Sec. II, the experimentally measured (and thus effective) resistivities, which relate the actual externally 
applied current density J to the electric field induced by vortex motion, considered in terms of its components along 
the directions both parallel and normal to the current. Here this model is reviewed and cast in tensor form, showing 
that it provides the same expression as Eq. (25) obtained from the force Equation. 

By isolating J T in Eq. (6), substituting it in Eq. (3) and considering that J s = J,sB, the following identity holds: 

pfjT=p { ff ) J-p ( ff ) BJ s (Bl) 
and by taking the dot product by B, considering Eq. (5), one obtains: 

Js = ) "J, ~ B (B2) 



■ B 



It is worth noting that J s and Jt are different from the components of J respectively parallel (J//) and perpendicular 
(J_i_) to B: in general Jt has itself a component parallel to B. Only when the superconductor is isotropic one has 
J s = J a and Jt = J±, as it can be checked considering Eqs. (5) and (B2). 
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Substituting Eq. (B2) together with Eq. (6) into Eq. (3) yields, with some tensor algebra: 

1 



E = 



(pfB) ■ B 



{Bx[($f)x($2)]} 



This equation allows to determine pg defined in Eq. (4). Since pg) is diagonal the above Equation can be rewritten 
by exploiting property (Al): 

)\ 

By comparing Eq. (4) and Eq. (B3), one can write down explicitly pg: 

p ff (6,<jy) = ~B x \ , -, .: . i /■■ 





=(») 

pff 






■ B 



=(») 
Pff 


=(0-i 
Pff 


(4°*) 


■B 



e 2 (M) 



(0 

%,11 



M 



■M" 1 IB, 



MB B 



(0 



M" 1 B* 



--M(0, 



(B3) 



(B4) 



which is equal to Eq. (25). 



Appendix C: The out-of-plane intrinsic flux flow resistivity cannot be measured with a single setup 

Without loss of generality I consider the static field B oriented according to the angles (9, 4>) and a current density 
J = J J E x-z plane, with unit vector J = (sin a, 0, cos a) (a is the angle between J and the z axis). Despite the 
choice J e x-z plane, the relative orientation of B and J is completely free, as long as a and <f> are free angles. In this 
configuration, the measured flux flow resistivity along J is, according to Eq. (26): 

J (a a \ Pff!nW 
p J (0,(j),a)- 



^(6)^6 + ^(6)^6 
(pg\i(9) sin 2 0sin 2 9 sin 2 a + p'g 33 (9)(cos 2 9 sin 2 a+ 
— 2 cos (j> cos 9 sin 9 cos a sin a + sin 2 9 cos 2 a)) = 

= pgjnW 

p^^smH + pf^cos^' 
[a 1 (<f>,6,a)p^ n (6)+a 3 (ct ) ,6,a)pg}] 33 (6)} (CI) 

where ai((j>, 6, a) and a 3 ((p,6,a) represent in compact way the combination of trigonometric functions appearing in 
the second equality. In order to obtain pg\ 3 (6) from the measured p J (9,(j),a) without being forced to know also 
Pfj\i(9)-> in the last equality p < g 11 (9) must disappear for each 9 with the freedom to choose properly <j> and a. This 
is impossible since the term in square parentheses should be equal to a rational function of pg n and Pg 33 - 
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On the other hand, it can be easily checked that it is possible to obtain pg ' n only: it is sufficient to choose <j> and 
a (e.g. maximum Lorentz force configuration: <j> — and a — ir/2) so that the term in square parentheses simplifies 
with the denominator of the fraction. 
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